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While the impurity-induced nanoscale electronic disorder has been extensively reported in the 
underdoped iron pnictides, its microscopic origins remain elusive. Recent scanning tunneling mi- 
croscopy (STM) measurements reveal a dimer-type resonant structure induced by cobalt doping. 
These dimers are randomly distributed but uniformly aligned with the antiferromagnetic a axis. 
A theory of the impurity-induced quasiparticle interference patterns is presented that shows the 
local density of states developing an oscillatory pattern characterized by both geometry and orbital 
content of the reconstructed Fermi pockets, occasioned by the pocket density-wave (PoDW) order 
along the b axis. This pattern breaks the d symmetry and its size and orientation compare well 
with the dimer resonances found in the STM experiments, hinting at the presence of a "hidden" 
PoDW order. More broadly, our theory spotlights such nanoscale structures as a useful diagnostic 
tool for various forms of order in iron pnictides. 



The iron pnictide high-T c superconductors [TJ [2] ex- 
hibit several remarkable features [3J. Among them 
are the proximity and interplay of antiferromagnetism 
(AFM) and a structural transition [H [5] , manifested by 
anisotropy in electronic properties [BTfTT] . Various exper- 
iments reveal strong correlation between these two tran- 
sitions [T2J[T3]. A prevalent explanation is that the struc- 
tural transition results from the fluctuations of incipient 
AFM order [TH [TS] . Alternatives include the key role of 
orbital degrees of freedom |16j and, in particular, the pro- 
posal that the structural transition originates from the 
pocket density-wave (PoDW) [T7J in parent compounds, 
a "hidden" order responsible for orbital ferromagnetism. 

Recent scanning tunneling microscopy (STM) exper- 
iments observe the nematic-type electronic formations 
developing around dopant atoms [9] HE] ■ These forma- 
tions appear as "dimer resonances," with two neighbor- 
ing peaks separated by ~ 6 — 8 lattice spacings. Im- 
portantly, the dimers are oriented along the a axis of the 
pnictides' unidirectional AFM order. In this Rapid Com- 
munication, we show that both the appearance of dimer 
resonances and their size and orientation provide a di- 
rect insight into the electronic structure and correlations 
in iron pnictides. In this regard, the properties of such 
dimers can serve as a diagnostic tool to unravel the na- 
ture of the underlying microscopic ground state and the 
sequence of Fermi-surface reconstructions left in the wake 
of various itinerant density- wave (DW) orderings [TU] . 

These are our main results: The dimer resonance [9lfl8] 
is a consequence of a type of "reconstructed" nesting [20] 
characterized by the wavevector q a s» (0.4-7T, 0). Such 
nesting tendency is manifest in our detailed calculations 
within the three orbital model [21], where q a - associ- 
ated with the short axis of the elliptical electron pocket 

emerges as a prominent feature of the reconstructed 
band structure, itself occasioned by the transition from 
the paramagnetic phase to the PoDW. Specifically: i) In 



the paramagnetic phase, the Fermi pocket is far from any 
nesting at q a . Consequently, no resonances appear; ii) As 
a PoDW is formed, the electron pocket is itself deformed 
while, simultaneously, sections of the reconstructed hole 
pockets become flatter as the PoDW order parameter in- 
creases. This gives rise to the "reconstructed" nesting 
at q a ; Hi) q a is perpendicular to the PoDW, producing 
the real-space dimer pattern breaking the C4 symmetry. 
The resonant pattern arises when randomly distributed 
dopant impurities induce two peaks in the local density 
of states (LDOS), both their separation and direction set 
by q a ; iv) When a PoDW and a SDW coexist, two nest- 
ing vectors q a ~ (0.47r,0) and q\ ps (0,0. Air) emerge from 
the reconstructed pockets. With e y and e x coupled to 
the inner and outer hole pockets, h% and hi, respectively, 
the susceptibility near q a remains much larger than its 
counterpart near q%. Again, the C4 symmetry is broken 
resulting in dimer patterns observed in [5] [TS] . 

We now supply the details behind the above physical 
picture. Experiments and theoretical calculations indi- 
cate that iron pnictides contain four disconnected Fermi- 
surface pockets (Fig. 1) [22] [23] • Additional parts of the 
Fermi surface are present in some materials [23] but this 
is not important for the physics discussed here. When 
two hole pockets closely match two electron pockets, as 
is the case in many iron pnictide parent compounds [23] , 
the geometric nesting favors the formation of two DWs 
[17] : The PoDW partially gaps e y and one hole pocket 
and induces a structural transition, while e x and the re- 
maining hole pocket form the partially gapped SDW. The 
theory [17] naturally explains the proximity of the two 
transitions and accounts for the observed orbital ferro- 
magnetism. 

While geometric nesting tendencies are important [22] , 
they are not the complete story: The orbital content of 
the Fermi pockets must be considered as well [24] . Here, 
we start with a simplified geometric model allowing for an 
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analytic glance at the physics and then fortify the results 
within a detailed numerical follow-up employing three- 
orbital content [21 j . known to capture the main features 
of real materials [53]. We model the Fermi pockets as: 



peak, implying no resonance around q a in the normal 
state of such a simple geometric model. These features 
of geometric nesting are echoed in a realistic calculation 
with full orbital content HUES]: 
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with m a > m > mj,, so that the two DWs are partially 
gapped when order parameters Ap Q rjw and A$dw are 
small. e y and e x bands are related by 90° rotation. To 
simplify the analytic calculation, we assume that hi and 
/i2 are isotropic and have the same dispersion relations. 




FIG. 1. (Color online) Fermi pockets in the unfolded Bril- 
louin zone of iron pnictides. Two hole pockets hi and hi are 
centered at the F = (0, 0) point. The elliptical electron pock- 
ets e x and e y are centered at M\ = (n,0) and Mi = (0, 7r), 
respectively, hi, hi, e x , and e y pockets exhibit strong nest- 
ing tendencies in many parent compounds. Here, we propose 
that the resonant electronic structure arises due to the "re- 
constructed" nesting with the wave vector ~ q a , the shorter 
axis of the e y pocket. 

It is known that static charge/spin susceptibilities con- 
tain peaks around the nesting vectors Mi = (vr, 0) and 
M2 = (0, 7r) in the above noninteracting model, implying 
a tendency to DW formation in moderately correlated 
iron pnictides [35]. Initially, by ignoring the orbital con- 
tent at the Fermi level within [25] , we can determine the 
analytic form for such geometric susceptibility at T = 0, 
within a single elliptical electron pocket ([!]): 
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x(<f) is effectively a constant for q/2 within the pocket, 
and decreases slowly as q/2 moves outside. Consequently, 
x(q) contains only an unremarkable ridge instead of a 
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where a, ft denote d orbitals and yn,s are the band in- 
dices. In the paramagnetic phase, the Green's func- 
tion (d^ 4 (Ki,a> n )dJ,.(fej-,a; n )) is finite only if ki = kj and 
fii = fij and thus Eq. (pi) gives 
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The factor of 2 in ^ is due to spin and n / (e) is the Fermi 
distribution function. Again, we find nesting peaks at 
(it, 0) and (0,7r) but no enhanced structure near q a . 

The situation changes once the DW order is estab- 
lished. Within a PoDW, which is the leading DW in- 
stability of Ref . [17] , the electron pocket e y couples with 
one of the hole pockets; in the following, we assume this is 
h%. The reconstructed Fermi surface is shown in Fig. [2] 
x(Qa) now increases, reflecting the nesting between the 
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FIG. 2. (Color online) As PoDW order develops, the Fermi 
surface reconstructs, and e y and the inner hole pocket hi 
are deformed. The curvature along e y decreases, leading to 
the "reconstructed" nesting at q a , the short axis of the e y 
pocket. This reconstructed nesting is absent in the paramag- 
netic phase, but becomes pronounced as A PoDW increases. 

reconstructed electron pockets: e y deformation promotes 
q a to a "reconstructed" nesting vector, connecting sepa- 
rate portions of two small e pockets, as depicted in Fig. 
[2] When Ap r)w is small, q = q a is a local maximum of 
X(o), with a strong peak when the nesting is optimized, 
i.e. (d 2 \/dk 2 y ) k ^ ±qj2 ky=Q = 0, A being the energy of 

the reconstructed h pocket. Within our simple model ([!]) 



- m b 



(4) 



3 



optimizes nesting. For Ap a p,w A opt , % increases by 
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Afc is width of the reconstructed h pocket, vp is the Fermi 
velocity || x, k = q a /2, and (3 = {d 4 e/dkf / ) k ^ =qj2k ^ Q . 




FIG. 3. (Color online) Left: Feynman diagrams for \ c in the 
PoDW phase. All other diagrams cancel out. Right: X c (k) 
with Ap dw ~ 20meV in the realistic three-orbital model. 
The original peaks at (0, ±7r) in the paramagnetic phase are 
suppressed by finite Ap pw- Instead, x(k) contains a plateau 
centered at q a m (±0.4-7r,0), breaking the C4 symmetry. This 
plateau arises from the "reconstructed" nesting at q a . 

Again, we fortify the above picture with the real- 
istic calculation of charge susceptibility x c {<i) within 
the three orbital model. In the PoDW state, 
(d Mi (fci, (kj,uj n )) is finite only if /x, = [ij and 

fcj = kj or (fj,i,Hj) = (e y ,h 2 ) and k l = kj + (0, tt). Fig- 
ure [3] displays Feynman diagrams contributing to x c (^)- 
Additional diagrams, composed of one normal and one 
anomalous Green's function, break time reversal symme- 
try and cancel out upon summation over fc. \ c IS ca l~ 
culated within the three-orbital model and is shown in 
Fig. [3] Evidently, a high plateau around (±0.47r, 0) arises 
due to the "reconstructed nesting" at vector q a . This is 
qualitatively consistent with the results in Eq. ([6| for 
the simple model (JlJ. In addition, the original peaks at 



M% = (0, ±7r) are suppressed due to the formation of a 
PoDW. 

At low T, Ref. [17] predicts the coexistence of a PoDW 
and a SDW, and therefore, leads to two "reconstructed" 
nesting vectors, q a and q bl related by 90° rotations. If 
Ap DW 7^ A-sdw, the C4 symmetry is already broken by 
this unequal pairing. In Ref. [T7| , A p DW > A5 p>w and 
this naturally leads to x(<la) > x{%)- When Ap p>w = 
Agp>w, however, the difference between x(<f a ) and xioV) 
could still arise from the distinct orbital content of Fermi 
pockets. To illustrate the effect of orbital content, we 
assume [53] 
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where 9 is the polar angle at the Fermi pocket. 
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FIG. 4. C4 symmetry is broken due to orbital content. Only 
the most significant diagrams contributing to the difference 
between x(<?f>) and x(9a) are displayed. The upper and lower 
two contribute to x(Qa) and x(<lb), respectively. This differ- 
ence arises from the vertex, which includes orbital overlaps be- 
tween incoming and outgoing fermions. The diagrams which 
have C4 symmetric partners are not shown. 

Figure. [4] shows the Feynman diagrams contributing to 
the susceptibility difference 5\ = x(Qa+M 2 )— x(qb+Mi). 
Quantitatively, S\ is dominated by the first diagram: 
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The vertex factors cos 2 9(k + q) and sin 2 9(k + q) come 
from the orbital component of the charge density p(r) 
when calculating x(<f a ) and x(Qb), respectively, 
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n/ is the Fermi function, and A is the energy of the recon- 
structed, partially gapped fermions. The other diagrams 
give similar contributions, with different vertices. 

It is now difficult to get an analytic expression for 
Sx even within the model ([!]). However, for Ap p>w — 
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with 5k, vj and (3 from Eqs. ([7|-(|9|. The anisotropy is 
determined by the choice of coupling between the e and 
h pockets: x(g a + M 2 ) > {<)x{9b + Ml) when e y and e x 
couple to h 2 (hi) and hi(h,2), respectively [55] . 




FIG. 5. (Color online) X c (k) with Ap Diy = 20meV and 
Asdw = 15meV within the three-orbital model. The original 
peaks around (0, ±7r) and (±7r, 0) in the paramagnetic phase 
are suppressed due to finite Ap pw- In addition, two small 
peaks emerge around (±0.47r, ±7r), breaking the C4 symmetry. 

Again, wc follow up with the realistic calculation of 
X C (<D based on the three orbital model. The results are 
displayed in Fig. [5] The peaks at the nesting vectors 
Mi and M 2 are suppressed due to the PoDW and SDW. 
Evidently, x(fc) breaks the C4 symmetry, as expected. 

STM experiments reveal an anisotropic electronic 
dimer structure developing when the iron pnictide is 
doped with cobalt [H] . Our theory of reconstructed nest- 
ing can account for this phenomenon, by considering the 
impurity-induced quasiparticle interference in the parent 
compounds. First, the local potential of the cobalt atom 
is H imp = Ea, a (V s a + aV£)dUr)d a (r), where V s and 
V m are the nonmagnetic and the magnetic parts of the 
impurity potential, respectively, and are given for each 
orbital a in [24]. The LDOS with a single impurity is: 
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Go is the bare Green's function matrix in orbital space. 

Figure. [6] shows the LDOS on a 10 x 10 square lat- 
tice, centered around the cobalt impurity. Note the 



FIG. 6. (Color online) The LDOS integrated from E = to 
37meV [l8| . with Ap dw = 20meV, Asdw = 15meV, and 
the impurity located at f = (0, 0). We assume e y (e x ) couples 
to h,2 (hi). LDOS peaks at ~ r — (±3a, 0), fixing the dimer 
size to 6a. The dimers are aligned with the x axis, reflecting 
the orientation of the peaks of charge susceptibility (Fig. [5]l . 
This anisotropy is due to different magnitudes of DW order 
parameters and different orbital content of individual pockets. 



pronounced peaks at r — (±3a,0), giving rise to the 
anisotropic dimer structure, observed in [18] . The 
dimer's size is tied to the magnitude of the "recon- 
structed" nesting vector q a « (0.47r,0). Since the cobalt 
impurity potential is repulsive, LDOS is small at the 
impurity site and large when r — (±ai27r/0.47r, 0) « 
(±3a, 0). Therefore, the size of the dimer is 6a, close to 
8a of Ref . [18 . In addition, the C4 symmetry is broken 
— even if Ap Dw = Asdw — due to the different orbital 
components within each pocket. The charge susceptibil- 
ity, therefore, peaks along the direction perpendicular to 
PoDW, and results in the a-axis orientation of the dimer. 

In summary, we have shown that the anisotropic elec- 
tronic dimer structure can be understood based on the 
"hidden" PoDW order in parent compounds. This order 
induces a "reconstructed" nesting vector q a , at which the 
charge susceptibility x(<f Q ) develops pronounced peaks. 
Furthermore, the C4 symmetry is genetically broken due 
to the orbital components. The dimer, therefore, points 
along the a axis of AFM order, in accordance with ex- 
periments. More generally, our results accentuate the 
potential of local probes as diagnostic tools in unraveling 
the patterns of various forms of order in iron pnictides. 
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